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We employ the Dynamical Cluster Approximation (DCA) 
in conjunction with the Fluctuation Exchange Approximation 
(FLEX) to study the Hubbard model. The DCA is a tech- 
nique to systematically restore the momentum conservation at 
the internal vertices of Feynman diagrams relinquished in the 
Dynamical Mean Field Approximation (DMFA). FLEX is a 
perturbative diagrammatic approach in which classes of Feyn- 
man diagrams are summed over analytically using geometric 
series. The FLEX is used as a tool to investigate the com- 
plementarity of the DCA and the finite size lattice technique 
with periodic boundary conditions by comparing their results 
for the Hubbard model. We also study the microscopic theory 
underlying the DCA in terms of compact (skeletal) and non- 
compact diagrammatic contributions to the thermodynamic 
potential independent of a specific model. The significant 
advantages of the DCA implementation in momentum space 
suggests the development of the same formalism for the fre- 
quency space. However, we show that such a formalism for the 
Matsubara frequencies at finite temperatures leads to acausal 
results and is not viable. However, a real frequency approach 
is shown to be feasible. 



I. INTRODUCTION 

Non-local correlations play an important role in the 
physics of strongly correlated electron systems such as 
high-Tc superconductors, heavy fermion metals, etc. The 
Dynamical Mean Field Approximation (DMFA) [1,2], in 
which all the non-local correlations are ignored, can cap- 
ture some of the major features of strongly correlated 
systems. Nevertheless, the non-local correlations become 
crucial in the physics of phases with non-local order pa- 
rameters such as d-wave superconductivity. Even phases 
with local order parameters such as commensurate mag- 
netism can be significantly affected by the non-local cor- 
relations (e.g. spin waves) ignored in the DMFA. 

The early attempts to extend the DMFA by including 
non-local correlations resulted in the violation of causal- 
ity which is a requirement for positive definiteness of 
the spectral weight and the density of states (DOS) [3]. 
The Dynamical Cluster Approximation (DCA) is a fully 
causal technique used to systematically add nonlocal cor- 
rections to the DMFA by mapping the lattice onto a 
self-consistently embedded cluster problem. The map- 
ping from the lattice to the cluster is accompanied by 
coarse-graining the lattice problem in its reciprocal space. 
Thus far, the DCA has been combined with Quantum 
Monte Carlo (QMC) [4,5], the Non-Crossing Approxima- 
tion (NCA) [6] and the Fluctuation Exchange Approxi- 
mation (FLEX) [11,12] to solve the corresponding cluster 



problems. 

The FLEX is a perturbative diagrammatic approach 
in which classes of Feynman diagrams are summed to 
all orders using geometric series. [7,8] Others [9,10] have 
employed the FLEX for finite size lattices with periodic 
boundary conditions. Due to the absence of contributions 
from some relevant diagrams, the FLEX is not capable 
of addressing the Hubbard model physics in the strong 
regime precisely. However, the main objective of this 
work is to make a comparison between the DCA-FLEX 
combination results and previous finite size lattice FLEX 
calculations. It is hoped that this study will lead to a 
better understanding of the DCA. 

We earlier [11] suggested a prescription to correctly 
implement the DCA technique in the thermodynamic 
potential. This prescription will be discussed from a 
different point of view using a more general argument. 
Based upon the Green function's exponential fall-off as 
a function of distance, we conclude that compact di- 
agrams (two-particle irreducible in the thermodynamic 
potential) are better approximated using the DCA than 
non-compact (two-particle reducible) ones. Hence, the 
DCA is applied to the compact diagrams only and non- 
compact ones are calculated explicitly using dressed non 
coarse-grained Green functions. 

In this work we also consider the extension of the DCA 
to frequency space. The many-body theory at finite 
temperatures is conventionally derived in terms of dis- 
crete imaginary Matsubara frequencies. We illustrate 
that even for a self-consistent algorithm like the FLEX, 
coarse-graining the imaginary frequency propagators re- 
sults in causality violations and can not be implemented. 
However, a real frequency formalism is shown to be causal 
and applicable not only to the FLEX, but also to other 
cluster solving methods such as the NCA. 

This paper is structured as follows. In the next three 
sections, we briefly review the DCA, its application to 
the Hubbard model, the FLEX, and then we describe 
how the FLEX and the DCA may be merged into a sin- 
gle algorithm. In the next three sections, we use the 
FLEX-DCA, in comparison to the FLEX for finite-sized 
systems, to explore the properties of the DCA. The last 
two sections, are devoted to a microscopic derivation of 
the DCA, and to an extension of the DCA to frequency 
space. 
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II. DYNAMICAL CLUSTER APPROXIMATION 
(DCA) 



which indicates that the momentum is conserved modulo 
Afc for transfers between the cells. 



Both the DCA and the DMFA may be derived by ex- 
ploring the momentum conservation in the diagrammat- 
ics. As depicted in Fig. 1, momentum conservation at 
each vertex is described by the Laue function: 



ki,q+k2 



(1) 



In the DMFA, momentum conservation at the internal 
vertices of irreducible Feynman diagrams is completely 
relinquished. I.e., the DMFA simply sets A = 1. [13] 
Hence, we may sum freely over all the internal momenta 
entering and leaving each vertex. Only local contribu- 
tions survive the sum. Thus, this is equivalent to map- 
ping the lattice problem onto a self-consistently embed- 
ded impurity problem. The DMFA becomes exact at 
infinite dimensions. [14] 

The DCA is an approach to systematically restore the 
momentum conservation relinquished in the DMFA. In 
the DCA, the first Brillouin zone in the reciprocal space 
is divided into Nc equal cells of linear size Afc labeled by 
K in their centers, and the momenta within each cell are 
labeled by k. Then k = K + k (c.f. Fig. 2). To visuahze 
this scheme in the real lattice, one could consider tiling 
the lattice of N sites by N/Nc clusters each composed of 
Nc = sites where L is the linear size of the subcell 
and D is dimensionality (c.f. Fig. 2 for L = 2). We 
will use this picture in section VIII while discussing the 
microscopic theory of the DCA. We label the origin of 
the clusters by x and the Nc intercluster sites by X. So 
for each site in the original lattice x = X -I- x. Care must 
be taken when choosing the cluster geometries in order 
to preserve the lattice point group symmetry and also 
satisfy some other criteria for cubic or square lattices. 
[15] 

In the DCA, we first make the following separation in 
Eq. 1 

A^Y^ g»(i+X)-(Ki-Q-K2+ki-q-k2) _ 

X 

The products Ki • x, Q • x and K2 ■ x = 2mr where n 
is an integer. Therefore, their associated phases may be 
neglected and Eq. 1 splits into two parts 



N 

'^kl,q-|-k2 ^'^ '5Ki,Q+K2 



(3) 



The DCA also ignores the phases e due to the posi- 
tion of the cluster in the original lattice and (far less im- 
portant) e"''' "^ corresponding to the position within the 
cluster. As a result, it approximates N/Nc S^^ ^_|_^^ = 1, 
so that 




FIG. 1. A typical vertex in a Feynman diagram with solid 
lines as one particle Green functions and the wiggly line as 
an interaction. 

The approximation made through the substitution 
A ^DCA corresponds to replacing all internal legs 
in the compact (skeletal) diagrams by the coarse-grained 
Green function G and interaction potential V defined by 



N 



G(K,c^„) = ^^G(K + k,c^„), 



N 



(5) 



(6) 



In section VIII, we will define the compact and non com- 
pact diagrams and elaborately discuss why only the com- 
pact ones undergo the coarse-graining approximation. 

Replacing A by Ajjca tremendously reduces the com- 
plexity of the problem because instead of having to per- 
form sums over all N states in the entire first Brillouin 
zone, we have sums over only a set of Nc states where 
Nc « N. 
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FIG. 2. The real lattice clusters (right) and (left) the first 
Brillouin zone divided into subcells. 



III. HUBBARD MODEL 

We will apply the DCA to study the Hubbard model 
Hamiltonian incorporating interactions between the elec- 
trons on a lattice. It includes a tight-binding part due 
to the hopping of electrons among the sites and an inter- 
action between the electrons. The general Hamiltonian 
reads 



2 



where 



and 



H = Ho + Hj, (7) 



CT <ij> 



iJ/ = i ^ njania'V{Rj - Ri) 

jLcrcr' 

V{R)^^. (9) 

Factors t and correspond to electron hoppings and 

Coulomb interactions respectively. Later in the paper, we 
will study the simplest Hubbard interaction which is fully 
local and only between electrons sitting at the same site 
having opposite spin directions. The interaction strength 
is a constant called U. Hence, for the local model, Eq. 9 
simplifies to 



(10) 



In terms of the vertex properties addressed in section H, 
since the interaction is local and therefore independent 
of q, we may sum freely over the q momentum for a pair 
of Laue functions in Eq. 1 sharing a common interac- 
tion wiggly line as depicted in Fig. 1. As a result, the 
corresponding Laue function will become 



N 



,ix-(ki-q-k2) 



E' 



,iy-(k3+q-k4 



A^^f^ki.q+ks (^ka+q.ki = N5\ 



ki+k3,k2+k4 



(11) 



and analogously for the DCA, by summing freely over Q 
A_DCA = Nc (5ki+K3,K2+K4 • (12) 



IV. FLUCTUATION EXCHANGE 
APPROXIMATION (FLEX) 

In the Feynman diagrammatics of the Hubbard model 
with a local interaction, all the interactions (wiggly lines 
in Fig. 3) contribute a c-number U from Eq. 10. The 
electronic Green functions (solid lines) which interact 
with one another should have opposite spins. Consid- 
ering these restrictions, had we been able to include all 
the possible diagrams in our expansion we would have 
solved the problem exactly. However, in practice this is 
not feasible. 

The Fluctuation Exchange Approximation (FLEX) 
was introduced as an approximate technique to simplify 
this diagrammatic sum [7,8] , while retaining a conserving 



approximation. In the FLEX, the interaction part of the 
Hubbard model Hamiltonian is treated perturbatively by 
selecting a certain class of all the possible diagrams which 
may be summed as a geometric series. Following Baym 
[16], we define the generating functional $[G(A;, tJn)] as 
the collection of all the selected families of diagrams il- 
lustrated in Fig. 3. Therefore, $[G(A:,w„)] for the FLEX 
can be written 



df 

ph 



-^Tr[xph]' + ^Tr[ln(l + Xph) 



Xph + -jX ph\ 



ph 



3 1 

- Tr [ln(l - Xph) + Xph + 2^%^ , 



1 



$pp = Tr [ln(l -f Xpp) - Xpp + pp] , 



(13) 

(14) 
(15) 
(16) 



where Tr = {T/N) with T the temperature and 

N the number of lattice sites. The particle-hole and 
particle-particle susceptibility bubbles are 

Xpp{q,U}n) = 

U{TIN) + + ^rn)Gi-k, -U^) , (17) 

k m 

Xph{q,OJn) = 

-U{T/N) 5] E + + ^™)G(fc, ^m) , (18) 

k m 



(a) t 



(b) t 



(c) 




(d) 



FIG. 3. (a) Lowest order diagram (Hartree term), (b) lon- 
gitudinal spin and density fluctuation $ diagrams with an 
even number of bubbles ($p{,), (c) transverse spin fluctua- 
tions ($p{) and (d) particle-particle fluctuations ("l>pp). 
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The self-energy and the Green function are defined by 
1 S^G] 



E(fc,w„) 



2 5G(fc,c^„) 



(19) 



(20) 



where G'"' is the non-interacting one particle Green func- 
tion defined by 



G(")(k,c^„) 



(21) 



with the non-interacting Hubbard model dispersion 
and /i the chemical potential. 

Calculating the self-energy for Eq. 13 using Eq. 19 we 
get 

I](fc,a;„) = UiT/N)J2Y. [V^^^Kl^^ra) x 



G{k - g, w„ - w„) - V'-PP^ {q, tJ„) X 
G(-fc + g, -cj„ + w„)] , 



(22) 



in which 



1 3 

fl + X..(.,-„)'^] + 2^^''('^'"'"»" 

1]- 



"1 - Xph{q,i^7n) 



Xpp{q,^7i 



-1 



(23) 



(24) 



Eq. 23,24 for the potential functions V^p'^^ and ^(pp) 
are geometric series for Xph and Xpp similar to the ran- 
dom phase approximation (RPA) results. The Hartree 
term contribution to the self-energy has not explicitly 
appeared in Eq. 22 as it is constant and can be always 
embedded in the chemical potential in Eq. 21. 

The difference A57(T, ^) between interacting and non- 
interacting thermodynamic potential functional is also 
expressible in terms of the Green functions, self-energy, 
and $[G] 



and to study finite sized systems with periodic boundary 
conditions, we can study the differences between these 
approaches. For example, as we will shortly illustrate, 
the complementarity of the DCA and finite size lattice 
techniques is manifest in the FLEX. The FLEX can also 
be invoked as a good test for the microscopic theory of 
the DCA and the coarse-graining effects in the compact 
and non-compact diagrams for Eq. 25. 



V. THE COMBINATION OF THE FLEX AND 
DCA (ALGORITHM) 

In the combination of the FLEX and DCA, our goal is 
to calculate the self-energy in Eq. 19 whereby we con- 
struct the dressed Green function for the lattice as a 
building block for all the relevant physical quantities. 
We start out with the bare (non-interacting) Green func- 
tion G^"-* (k, z) defined in Eq. 21 with z the Matsubara 
frequency (complex). We coarse-grain G as directed in 
Eq. 5 and calculate the self-energy using Eq. 19. This is 
used to recompute the dressed Green function 



G(k,z) 



1 



Ek + /-i - Sdca(K, z) 



(26) 



where the index DCA in E£)ca(K, z) indicates that we 
have used coarse-grained G for the construction of self- 
energy. The new G is coarse-grained and used to calcu- 
late a new estimate of T,dca{^, z). We repeat this pro- 
cess iteratively until convergence at a desired tolerance is 
obtained. The final self-energy is used to construct the 
dressed Green function in Eq. 26, required to compute the 
physical quantities such as spectral function, the density 
of states (DOS), etc. The algorithm of this calculation is 
demonstrated in Fig 4. 



Solve Cluster 
Problem 

(FLEX) 



Z(K) 



Nc 

G(K) = ^Z 



1 



An{T, n)^n-no^ -2Ti[J:G - ln(G/Go)] + $[G] . 

(25) 

In the FLEX, since we include only a limited set of all 
the diagrammatic contributions, we do not anticipate to 
precisely address the Hubbard model physics. However, 
there are a number of significant physical features such as 
anti-ferromagnetic order at half filling and low tempera- 
tures that this approximation is able to capture. More- 
over, by using the FLEX both together with the DCA 



FIG. 4. The FLEX-DCA numerical algorithm. The dressed 
G is solved self-consistently with E. The iteration process 
stops whenever convergence is achieved. 



VI. COMPLEMENTARITY OF THE DCA TO 
THE FINITE SIZE LATTICE APPROXIMATION 
WITH PERIODIC BOUNDARY CONDITIONS 

In the half-filled Hubbard model, the antiferromagnetic 
correlation length ^ increases with decreasing tempera- 
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ture and diverges at the phase transition. In a finite size 
lattice with periodic boundary conditions, as the tem- 
perature drops, once the correlation length reaches the 
size of the lattice, the system is fully frozen and there is 
a gap to excitations (c.f. Fig 5. a). In contrast, in the 
DCA, the correlations are confined within clusters of size 
Nc << N (the size of the entire lattice) and they never 
reach the size of the lattice. As we lower the temper- 
ature, the correlation length approaches the size of the 
cluster but since the lattice remains in the thermody- 
namic limit, it never freezes (c.f. Fig 5.b). By increasing 
the size of clusters in the DCA, we take longer ranged 
correlations into account so the gap will become more 
pronounced. Consequently, correlation induced gaps are 
generally overestimated in the finite size lattice, while in 
the DCA they are underestimated. 




• f • 




Cluster 



(a) 



(b) 



freezing 



no freezing 



FIG. 5. (a) The finite size lattice with periodic boundary 
conditions, size L and correlation length ^ and (b) the lattice 
with clusters of size L and the correlation length ^. 
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FIG. 6. The Green function G{k = (0,7r),r) at T = 0.10 
and U/t = 1.57 versus imaginary time t/P {13 = l/T). Finite 
size results (filled symbols) and DCA results (open symbols) 
are displayed for lattice sizes of A'^ = 8 x 8, 12 x 12, and 
20 X 20. The curves show the complementary approach of the 
two methods to the thermodynamic limit. 

This complementary behavior may be seen in G'(fc,r), 
with k = (0, tt) computed using finite size lattices with 
periodic boundary conditions and the DCA. As illus- 
trated in Fig. 6, by increasing the size of the finite size 
lattice and the DCA cluster, the Green functions con- 
verge from opposite directions. In the finite size lattice, 
the Green function (e.g. at r = /3/2) decreases with the 
increase of size which is consistent with overestimating 
the gap; while in the DCA, the Green function increases 
as the cluster size grows consistent with underestimating 
the gap. It is also observed that the convergence in the 
DCA is much faster, meaning that the result of the DCA 
is closer to the true curve at a given cluster size. Both 
finite size lattices with periodic boundary conditions and 
the DCA converge with corrections of 0{X/L^) with L 
the linear size of the finite lattice or the DCA cluster 
and A a coefficient. [17] The faster convergence of the 
DCA corresponds to its smaller A compared to finite size 
lattices. 
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FIG. 7. The Green function G{k = (0,7r),r) at r = (3/2 
(left) and t = /3/4) (right) versus 1/L^ for the DCA and 
finite size results in Fig. 6. The linear extrapolations meet at 
a point representing the thermodynamic limit of the Green 
function of the Hubbard model evaluated within FLEX. 

In Fig. 7, the values of the G{k = (0,7r),T = (3/2) 
and G{k = (0, 7r),T = /3/4 in Fig. 6 have been plotted 
versus 1/L^ for both the DCA and finite size results. The 
Green function behaves linearly as a function of 1/ for 
large L. The extrapolations of the DCA and finite size 
results meet as L — > cxd, approximating the value of the 
dressed Green function in the thermodynamic limit. The 
complementarity of DCA and finite size methods allows a 
determination of the thermodynamic limit of imaginary 
time Green functions with unprecedented accuracy. 
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VII. FINITE SIZE VERSUS THE DCA FLEX 
RESULTS FOR THE TWO DIMENSIONAL 
HUBBARD MODEL AT HALF-FILLING 

The Hubbard model at half filling undergoes a phase 
transition to anti-ferromagnetic order at low tempera- 
tures. According to the Mermin-Wagner-Hohenberg the- 
orem, for dimension D = 2 the critical temperature is 
zero. However, as we continue to lower the tempera- 
ture, close enough to zero, a pseudogap will appear in 
the density of states (DOS) as a precursor to the anti- 
ferromagnetic phase (which has a full gap as its signa- 
ture). An approach towards non- Fermi- liquid behavior 
is also visible in both the real and imaginary parts of the 
retarded self-energy. [9] 

In Fig. 8 and Fig. 9, the densities of states (DOS) for 
lattices with finite sizes of 32x32, 64x64 and periodic 
boundary conditions and coarse-graining cluster sizes of 
16x16 and 32x32 are plotted. We analytically continue 
the Green function in order to calculate the spectral func- 
tion A{k, ui) and DOS using the Fade approximation. [18] 
In this approximation, we generate a continued fraction 
interpolating all the data points and use it as an an- 
alytic function of the Matsubara frequencies a;„. The 
analytic continuation is accomplished by substituting a;„ 
with Lu + ir) where rj is a small positive shift. However, 
the errors inherent in the numerical Fourier transform 
(FFT) and also the sharp high-frequency behavior of 
the Green function, FLEX potentials and the self-energy, 
limit the accuracy of the Fade. The high-frequency be- 
havior is improved by implementing a more authentic 
cut off scheme introduced by Deisz et al. [19] in which 
the high-frequency tails of these quantities are Fourier q 
transformed analytically prior to any numerical FFT and '-' 
added back to the FFT outputs afterwards. In addition 
to this high-frequency cut off improvement, we also check 
for the analyticity of the Fade results in the upper-half 
frequency plane as a requirement for retarded physical 
quantities. This task is carried out by converting the 
continued fraction in the Fade into a ratio of two poly- 
nomials. The complex roots of these two polynomials are 
obtained via the Jenkins- Traub root finder routine. [20] 
Those orders of the Fade for which there exist poles in 
the upper-half plane are omitted unless these poles are w 
canceled by the roots of the numerator. The acceptable Q 
Fade results correspond to the highest order with no un- 
compensated poles in the upper-half plane. 

As seen in Fig. 8, by increasing the size of a finite lat- 
tice, the pseudogap occurs at higher temperatures and it 
also becomes less pronounced (sharper) as we approach 
the actual size of an infinite real lattice. The DCA yields 
a complementary behavior as shown in Fig. 9. By in- 
creasing the coarse-graining cluster size, similar to the fi- 
nite size lattices, the pseudogap is shifted towards higher 
temperatures. However, unlike the finite size lattices, 
for the DCA the precursor becomes more pronoimced 
(broader) as the cluster increases in size because the size 



of the lattice remains constant and the correlations are 
limited to the cluster size. Thus comparatively, the DCA 
underestimates the gap while the finite size calculation 
overestimates it. 

By comparing the results in Fig. 8 and Fig. 9 one may 
see that the 32 x 32 DCA cluster yields more realistic 
physics than the corresponding 32 x 32 finite size lattice. 
The 64 X 64 finite size lattice results are also close to 
those for the 32 x 32 DCA cluster at slightly lower T 
(eg, T = 0.055 for the finite size and T = 0.033 for the 
DCA). However since the sizes of clusters are consider- 
ably smaller than the sizes of lattices, the DCA signifi- 
cantly reduces the complexity of the problem and conse- 
quently the CPU time. In terms of the CPU time, the 
FLEX with the numerical Fourier transforms scales as 
AflnAf where Af is the product of the total number of 
Matsubara frequencies and the k points in the first Bril- 
louin zone. Hence, using a 32 x 32 cluster in place of 
a 64 X 64 lattice both having 1024 Matsubara frequency 
points roughly reduces the CPU time by a factor of 4.4. 
If the DCA cluster size Nc equals the size of the finite 
lattice N, the DCA requires somewhat more CPU time 
than the finite size lattice due to course-graining. Nev- 
ertheless, comparatively, for large finite lattices such as 
64x 64, the lattice size contributions to the CPU time sig- 
nificantly dominate the coarse-graining ones in a 32 x 32 
cluster. Thus, the 32 x 32 DCA cluster is much faster. 
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FIG. 8. The FLEX density of states (DOS) for a 32 x 32 
finite size lattice with periodic boundary conditions (top) ver- 
sus energy for U/t= 1.57 at various temperatures, and for a 
64 X 64 finite size lattice (bottom). The evolution of a week 
pseudogap is observed at T = 0.067 for the 64 x 64 lattice 
higher than T = 0.055 for the 32 x 32 one and the pseudogap 
is also much broader for the 32 x 32 lattice (128 x 128 finite 
size lattice results were obtained by Deisz et al. [9]). 

The FLEX often has difficulty reaching low tempera- 
tures. This is due to the fact that the Xph defined in 
Eq. 18 approaches unity as the temperature drops which 
in turn causes the V'^'^^^ in Eq. 23 to diverge. In the DCA, 
Xph approaches unity more slowly, allowing the calcula- 
tions to reach lower temperatures. One has to note that 
in the FLEX, the Xp^^ is defined as follows 



^^^(Q,o;„) = 



-U{T/N^) ^ ^ G(K + Q, a;„ + u;„)G(K, w„) , (27) 

K m 

with G defined in Eq. 5. Fig. 10 illustrates the saturation 

of Xph for both the DCA and finite size lattices. The Xph 
for the 32 x 32 finite size lattice (filled circles) saturates 
at higher temperatures compared to the 32 x 32 DCA 
cluster (open diamonds) indicating that in the DCA, for 
a certain cluster size, the precursor to the phase tran- 
sition can evolve to lower temperatures compared to a 
finite size lattice of the same size with periodic boundary 
conditions. However, for the correlation length £, > L 
the DCA approximation breaks down and replacing the 
self-energy by its coarse-grained counterpart is no longer 
accurate. Here, the DCA takes on significant mean-field 
character. 
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FIG. 9. The FLEX density of states (DOS) for a 16 x 16 
DCA cluster (top) versus energy for U jt = 1.57 at various 

temperatures, and for a 32 x 32 DCA cluster (bottom). For 
the 32 X 32 cluster, the evolution of a weak pseudogap is 
observed at T = 0.055 higher than T = 0.025 for the 16 x 16 
one and the pseudogap is also broader for the 32 x 32 cluster. 

Another feature of the Hubbard model near half-filling 
verified by the FLEX [9] is non-Fermi-liquid behavior. 
Here, this is studied by increasing the electron electron 
interaction U at a constant temperature. In Fig. 11 and 
Fig. 12 the real and imaginary parts of the self-energy 
at the X point (on the non-interacting Fermi surface) 
have been plotted versus energy for finite size lattices and 
the DCA respectively. As the interaction is increased, 
the negative slope in the real part turns positive around 
oj = which is inconsistent with the requirement that 
the renormalization factor [1 — 9i?eS(kir, a;)/9a;|aj=o]~^ 
should be smaller than unity in the Fermi-liquid theory. 
There also appears an anomalous inverted peak in the 
imaginary part at w = 0. 

As presented in Fig. 11, by increasing the length of the 
finite size lattice, the sharpness of the non-Fermi-liquid 
features is reduced. The same features for the DCA in 
Fig. 12 are slightly less pronounced and in a complemen- 
tary fashion to the finite size lattices, their sharpness is 
enhanced by increasing the size of the cluster. Thus, 
again the DCA underestimates the non-Fermi-liquid fea- 
tures while the finite size calculation overestimates it. 
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FIG. 10. The FLEX particle-hole bubble (xph) at 
Q = (tt, n) and Um = for a 16 x 16 and 32 x 32 DCA clusters 
(diamonds) and also for 32 x 32 and 64 x 64 finite size lattices 
with periodic boundary conditions (circles) versus tempera- 
ture T at U/t = 1.57. The relatively rapid saturation of the 
finite size lattice Xph compared to the more gradual approach 
of the DCA results towards unity is manifest (128 x 128 finite 
size lattice results were obtained by Deisz et al. [9]). 

Earlier in Fig. 10 it was shown that the FLEX particle- 
hole bubble Xph{q, t^m = 0) at g = Q = (tt, tt) approaches 
unity as the temperature is lowered. This causes the spin- 
fluctuation T matrix 



3r Xp/i(g,C^ 

2k- 



(28) 
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which is just the third term in y (p'') ((jf, lj^^) in Eq. 23 to 
peak around the [Q.ujm = 0) point. For real frequen- 
cies, Tu^a{Q,uj) has a sharp peak around u; = 0. Since 
Ta,a{Q,'^) is only used to construct the irreducible self- 
energy, within the DCA it is constructed from coarse- 
grained Green functions. Tliiis, the DCA counterpart of 
Eq. 28 is obtained by only replacing the Xph with Xpl^^ 
defined in Eq. 27. Fig. 13 and Fig. 14 show how this peak 
sharpens as the temperature decreases or interaction in- 
creases for finite size lattices and the DCA respectively. 
In the 64 x 64 finite size lattice (c.f. Fig 13 bottom), 
the peak continues to develop as the temperature is low- 
ered and the interaction is raised. At T = 0.067 and 
U = 1.6 where there exists a pseudogap in the DOS, 
the peak undergoes a significant growth compared to the 
other graphs shown in the same figure. The 32 x 32 
(c.f. Fig 13. top) lattice presents the same behavior with 
a slightly sharper but shorter peak. 

The DCA also illustrates the same type of peaks at a 
slightly higher interaction and lower temperatures for the 
16 X 16 cluster (c.f. Fig 14 top). Increasing the cluster 
size to 32 x 32 gives rise to higher peaks similar to the 
finite size lattices case but unlike the finite size lattices 
peaks become sharper as the size is increased. 
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FIG. 11. The FLEX real (top two) and imaginary (bot- 
tom two) parts of the self-energy at the X point for 32 x 32 
and 64 x 64 finite size lattices with periodic boundary con- 
ditions versus energy for T = 0.10 and two different interac- 
tions U . The inverted peak at = in the imaginary part 
and the positive slope in the real part are both signatures 
of non-Fermi-liquid behavior. For smaller lattice sizes these 
signatures axe more pronounced (128 x 128 finite size lattice 
results were obtained by Deisz et al. [9]). 



3 

w 
<u 
nc 



3 

w 
<u 
nc 



0.1 r 

0.05 


-0.05 
-0.1 

0.1 
0.05 


-0.05 
-0.1 



DCA 
Nc=16*ie 



T=0.10, U=1.0 
T=0.10, U=1.7 
T=0.10, U=2.1 





CO 





0.25 




0.15 


3 


0.05 


x" 


-0.05 


w 




E 


-0.15 




-0.25 




-0.35 



DCA 




- T=0.10, U=1.0 


Nc=32*32 




- T=0.10, U=1.7 






— T=0.10, U=2.1 








2 -1 





1 s 




0) 




DCA 




- T=0.10, U=1.0 


Nc=16*16 




- T=0.10, U=1.7 






- T=0.10, U=2.1 









-10 




CO 



10 



-10 




CO 



10 



8 





0.25 




0.15 




0.05 


t 


x" 


-0.05 


w 




E 


-0.15 




-0.25 




-0.35 



DCA 


T=0.10, 


U=10 


Nc=32*32 


T=0.10, 


U=1.7 




T=0.10, 


U=2.1 









-10 




CO 



10 



FIG. 12. The FLEX real (top two) and imaginary (bottom 

two) parts of the self-energy at the X point for 16 x 16 and 
32 X 32 DCA clusters versus energy for T — 0.10 and two 
different interactions U. The non-Fermi-liquid features occur 
at slightly higher interactions compared to the finite size lat- 
tices. Upon increasing the cluster size, these features become 
more pronounced (complementary to the finite size lattices 
with periodic boundary conditions). 

All the results illustrated in this section indicate the 
complementarity of the DCA to the finite size lattice 
scheme. It is also observed that the DCA is capable of 
reproducing relatively the same physics as the finite size 
FLEX at slightly different parameters but a lower CPU 
cost. The combination of these two facts makes this tech- 
nique a good candidate to be employed in the numerical 
treatment of a wide range of many-body problems. 
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FIG. 13. The FLEX real part of the T-matrix at Q = (tt, tt) 
versus energy for 32 x 32 (top) and 64 x 64 (bottom) finite size 
lattices with periodic boundary conditions. For the 32 x 32 
finite size lattice at T = 0.067 and U/t = 1.6 (solid), we 
observe a huge peak duo to the formation of a pseudogap in 
the DOS. By increasing the lattice size to 64 x 64 a higher 
and broader peak occurs at the same temperature. 
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FIG. 14. The FLEX real part of the T-matrix at Q = (tt, tt) 
versus energy for 16 x 16 (top) and 32 x 32 (bottom) DCA 
clusters. Peaks similar to the finite size lattices appear at 
slightly higher interactions and lower temperatures. By in- 
creasing the cluster size, the peaks grow higher (similar to 
finite size lattices) and also sharpen (complementary to the 
finite size lattices). 

Lastly, to better understand the effect of the DCA clus- 
ter embedded in a Fermionic bath, we can rewrite the 
course grained Green function defined in Eq. 5 as 



G{K,z) 



1 



z-eK-S(K,z)-r(K,z) ■ 



(29) 



where €k = Nc/N J2k ^K+k ^^'^ r(K, z) is the host func- 
tion. Maier et al. [6], define 



whereby r(K, z) can be expressed as 



N Z^k ' K+k 



G(K + k. z) 



i + ^EicWkGlK + k,^) 



(30) 



(31) 



By Taylor expanding t^-i-k ai'ound the cluster points 
K it is found that tj^^^ - 0{Ak) with Afc = 2tt/L. 
Thus, Eq. 31 yields r(K) ~ 0{{Ak)^) as Ak 0. To 
illustrate this, we calculate r(r = 0,r = 0) by sum- 
ming over all the K points and frequencies and plot 
it versus (Afc)^. Fig. 15 illustrates this linear behavior 
for Nc > 16. Nc = 1 holds complete mean field charac- 
ters and no non-local fluctuations. N^, = 4 is anomalous 
as explained in an article by Belts et al. [15] There, the 
finite size cubic lattices with less than six (four in 2D 
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square lattices) distinct nearest neighbors per each site 
are not used in finite size scalings for estimating the phys- 
ical properties of models like the spin one half XY ferro- 
magnet or the Heisenberg antiferromagnet. For = 4, 
because of the periodicity of the clusters, each cluster 
point is surrounded by two identical nearest neighbors 
in every direction and therefore has only two distinct 
nearest neighbors (c.f. Fig. 16). Thus, the effect of fluc- 
tuations are overestimated. [21] For Nc > 4 there is no 
such anomaly and hence, all the points present the linear 
behavior proven above. Nevertheless, calculations with 
A'c = 4 do a reasonable job in capturing the qualitative 
effects of corrections to the DMFA. 
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FIG. 15. The host function r(r = 0,t = 0) versus (AA:)^ 
The linear behavior beyond ATc = 4 is manifest, see also the 
inset. 
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FIG. 16. DCA clusters in the real space lattice with 
Nc = 4. Due to the periodicity of the clusters, points 2 and 4 
in squares and 1 and 3 in circles are equivalent and therefore 
point observes its nearest neighbors twice in each direction. 



VIII. MICROSCOPIC THEORY OF THE DCA 

In section IV we defined the thermodynamic poten- 
tial functional difference An{T, n) in terms of the Green 
function, self-energy, and $[G] (c.f. Eq. 25). In Eq. 13, 
includes all the compact (skeletal) Feynman dia- 
grams and the rest incorporates the entire non-compact 
contribution. [22] Typical compact and non-compact di- 
agrams are illustrated in Fig. 17. The non-compact dia- 
gram (a) consists of two self-energy pieces a and a' con- 



nected with two one-particle Green functions. Remov- 
ing these two Green functions would split the diagram 
into two separate pieces. In the compact diagram (b) 
two vertex parts F and F' with four Green functions are 
connected together. One can not split this type of dia- 
grams into two separate pieces by just removing two one- 
particle Green functions. As mentioned earlier, in the 
DCA, we employ coarse-grained Green functions to con- 
struct only the compact diagrams. The Green functions 
in non-compact diagrams are calculated directly using 
Eq. 26 in which the self-energy Sdca(K, z) is coarse- 
grained (the circles at the top and the bottom of the 
non-compact diagram in Fig. 17. (a)) 
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G(k) 
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G(k) G(k) 
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G(k+k-q) 



r'(k,k',q) 



FIG. 17. (a) typical non-compact (nan- skeletal) and (h) 
typical compact (skeletal) diagrams. 

Earlier in a shorter article for the Hubbard model in 
particular [11], we showed both analytically and numeri- 
cally that the error produced by coarse-graining the non- 
compact diagrams is significantly larger than the error 
produced by coarse-graining the compact ones. Here, we 
would like to give a more general argument in real space. 
We wish to emphasize two points in this new approach. 
First, since the derivation of the DCA in this section re- 
lies only upon the exponential fall off of the Green func- 
tion as a function of distance, it is far more intuitive than 
the momentum space argument in Ref. [11]. Second, it 
ties the derivation of the DCA to the original derivation 
of the DMFA in the limit of infinite dimensions, where 
similar arguments are employed. [14] 

The exponential fall off behavior occurs naturally in 
high dimensions. In the tight-binding Hamiltonian (c.f. 
Eq. 8), the factor t corresponds to the hopping of elec- 
trons among nearest neighboring sites. Thus, one could 
show that the real space Green function G{r) (we drop 
the frequency label from this point on for simplicity) for 
r nearest neighbor hops is proportional to G(r) ~ V as 
f — » 0. On the other hand, Metzner et al., and Miiller- 
Hartmann [14,13] have shown that in D dimensions, the 
factor t should be renormalized as tj^fD in order to have 
a finite density of states width as £) — > oo. As a result of 
this renormalization 

G(r) r^f r-. [l/^/DY ~ D-^'^ = e-''/'^^ 
meaning that G{r) falls off exponentially as a function of 
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In the DCA, we attempt to minimize the error due 
to coarse-graining the Green function (and potentials) 
in the Feynman diagrams. Consider the first non-trivial 
correction to the coarse-grained non-compact diagrams 
generated by replacing the explicit coarse-grained Green 
function lines by the non-coarse-grained ones as illus- 
trated in Fig. 18 



N, 



K 



a{K)a'{K)G{Kf 



(33) 



where K are the coarse-graining cells momenta and k = 
K -h k include all the momenta in the first Brillouin zone 
shown in Fig. 2. In this derivation we also presume that 
the self-energy is k independent and the entire k depen- 
dence is embedded in the Green functions. 





FIG. 18. First correction by non-compact diagrams, 



By breaking up the sums over k in Eq. 33 into K and 
k sums and writing all the G in terms of their Fourier 
transforms at the same time we get 

,5(i)[Af]„ep] ~ ^^a(K)a'(K) ^ G(xi)G(x2) x 



k,K 



g-ik.(x,+x.)g-iK.(xi+x.) _ iV ^ a(K)a'(K) 



K,ki,k2 



G(xi)G(x2)e-'''i-^ie-'''=-^=e-'^-(^i+^=) 



(34) 



in which we used Eq. 5 for G(K). According to Fig. 2, 
we can split x = X + x where x extend between two 
separate clusters while X always remain within a single 
cluster. By making this separation in Eq. 34 one picks 
up phases including products of — iK-X, — iK-x, — ik-X 
and — ik-x in their exponents. The phase associated with 
the product K • x = 2n7r with n an integer equals unity. 
The phases involving — ik- X products are also neglected 
as discussed in section II. Hence, Eq. 34 can be rewritten 
as follows 



a'{K) E E <^(^i+^i)^(^2+X2)e-^^-(^^+^^) X 

Xi,X2 Xi,X2 

E e-'^^ '^'e-'^^-''^ . (35) 

ki,k2 

Implementing the following substitutions 



E' 



-ik-(xi-|-X2) 



N. 



and E 

ki,k2 



iki-xig-ik2-X2 



= (-^)^<5ii,0<5i2,0 ■ 



(36) 



Eq. 35 simplifies into 

^W[Af]„ep]~^E^(K)^'(K)[ E G(Xi+x) 



K 



Xi,X2,X 

G(X2 - x)e-'K-(Xi+X2) _ J2 G(Xi)G(X2) x 



Xi ,X2 



iK-(Xi-|-X2)l 



(37) 



K 



G(X2+i2)e 



-jK-(Xi+X2 



'E' 



Xi,X2 Xl,X2 

ik-(xi-|-X2) 



Setting cr(K)cr'(K) = ^(K) and performing the K sum- 
mation 

5W[Al^„ep]~ E C(Xi+X2)[EG(Xi+x)x 

Xi,X2 5c 

G(X2-x)-G(Xi)G(X2)] = E C(Xi+X2)x 

Xi,X2 

E G(Xi + x)G(X2 - x) . (38) 

x#0 

Knowing that ^(Xi + Xa) = o-(X)cr'(X + Xi+ Xa) 
and also that the lowest order of cr(X) ~ G'^(X) ~ 
g-3|x|/rs conclude that in Eq. 38, the largest con- 
tribution is due to terms having Xi = — X2 or in other 
words, local ^. As shown in Fig. 19, the first term in the 
X sum corresponds to |x| = L (size of the cluster) and X 
can be as large as X = —(L — 1) in the opposite direction. 
Hence, the leading order term in Eq. 38 falls off as 

(^(1) [Aflncp] ~ C(0) X 2i:> X G(L ~ {L- 1)) X 

G(-L + {L- 1)) ~ 2D^{0)e^e^ = 215^(0)6^ , (39) 

where 2D is the number of |x| = L contributions in D 
different dimensions of a £> dimensional cubic lattice. In 

Eq. 39 we also used the fact that due to the lattice sym- 
metry, G(-X) = G(X). As £> 00, using r, = 2/lnD 

(5(i)[An„ep]~2£>^(0)e-'"^ = 2^(0)£'£>-^~ 0{l) . (40) 

which indicates the existence of non-local corrections to 
the non-compact contribution of the thermodynamic po- 
tential even at infinite dimensions. 
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FIG. 19. Two adjacent clusters with their X and x vectors. 

Now we replace the coarse-grained self-energy bubbles 
in Fig. 18 with coarse-grained vertices having four exter- 
nal legs and look at the difference between compact dia- 
grams with and without coarse-grained Green functions 
which are explicitly shown in the figure. Since we ear- 
lier dropped the frequency labels in the Green functions, 
here wc use indices 1, 2, 3 and 4 to emphasize that these 
Green functions have different frequency labels. The first 
correction to the compact contribution of the thermody- 
namic potential depicted in Fig. 20 is 

5(i)[Af2ep] ~ 1^ E r(Ki,K2,Q)r'(Ki,K2,Q)Gi(ki) x 



r'(Ki, K2, Q)e~'^^'(-^^+-^*)e~*^^'^"^^~'~"^*^e~**^'^"^^~"^*^ 

Er(x,x',x")x 

x,x' 
x" 

r'(x + Xi -t- X4, X' + X2 + X4, x" + X3 - X4) . 

Once again, the largest contribution is associated with 
local A, i.e., Xi = — X4, X2 = — X4 and X3 = X4. 

Therefore 

5(1) [A^cp] ~ E A(0) [ E Gi (-X - x)G2 (-X - x) X 

X X 

G3(X + x)G4(X + x) -Gi(X)G2(X)G3(X)G4(X)] = 
A(0) Gi(X + x)G2(X + x)G3(X + x)G4(X + x). 



(43) 



X.XT^O 



(44) 



kl,k2 



G2(k2)G3(q)G4(ki +k2 - q) - ^ r(Ki,K2,Q) x 



Considering |x| = L and X„ij„ = — (L — 1) 

S^^^A^cp] ~ A(0) x 2D x Gi(-L + {L - 1)) x 
G2{-L +{1- 1))G3(L - (L - l))G4(i: -{L- 1)) = 

^ ;) (45) 



2DA(0) 



2A(0)i:'"^ 



InD' 



r'(Ki,K2, Q)Gi(Ki)G2(K2)G3(Q)G4(Ki + K2 - Q) . 



(41) 



where similar to Eq. 33, all the vertices are coarse-grained 
but the Green functions are not. 
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FIG. 20. First correction by compact diagrams, ^''^^fAficp]. 

Following the same procedure as for the non-compact 

difference wc arrive at 

,5W[Afiep]~ A(Xi+X4,X2 + X4,X3-X4) X 

^1 1^2 

Gi(Xi - i)G2(X2 - i)G3(X3 + i)G4(X4 +S<i)- 

X 

Gi(Xi)G2(X2)G3(X3)G4(X4)] , (42) 
where 



which vanishes as D ^ 00. 

Comparing Eq. 45 with Eq. 40 shows that the first 
correction to the compact contribution of the thermo- 
dynamic potential falls off exponentially twice as fast as 
the equivalent correction in the non-compact contribu- 
tion. In addition, even at the infinite dimensional limit, 
there are corrections of order one to the non-compact 
contribution whereas for the compact diagrams the DCA 
becomes exact and there are no corrections. This justi- 
fies coarse-graining only in the compact diagrams. The 
Green functions in the non-compact diagrams have to be 
explicitly constructed from the coarse-grained self-energy 
E/3C'A(K,a;„) using Eq. 26. 
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FIG. 21. The compact contribution to the thermodynamic 
potential versus (Afc)^ at U/t = 1.57 and various tempera- 
tures, using coarse-grained Green functions. Tlic deviation 
from linearity at the lowest temperatures hint at the correla- 
tion length exceeding the cluster size. 
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FIG. 22. Non-compact contribution to the thermodynamic 
potential versus (Afc)'^ at U/t = 1.57 with and without 
coarse-grained Green functions. Using coarse-grained Green 
functions can result in a non-compact contribution with the 
wrong sign. 

In Fig. 21, the compact contribution of the thermody- 
namic potential difference constructed by coarse-grained 
Green functions is plotted versus (Afc)^ = l/L^. The 
variation of AOcp over the entire range is about 

1%. At very low temperatures (T < 0.066 in the inset), 
some deviation from linearity is observed due to the cor- 
relation length exceeding the cluster size L and therefore 
the approximation of E(A;) by T,£)ca{K) begins to break 
down. Fig. 22 illustrates the non-compact contribu- 
tion sketched versus l/L^ using both coarse-grained and 
non-coarse-grained Green functions. Using non-coarse- 
grained Green functions, the variation of AQncp over the 
entire range is roughly 30%. Coarse- graining the 

Green functions in these diagrams can even change the 
sign of this non-compact contributions, clearly indicating 
that coarse-graining the Green function is an unreahstic 
approximation for non-compact diagrams. However, one 
notices that for large cluster sizes the coarse-grained and 
non-coarse-grained results approach each other as the ap- 
proximation to the infinite lattice becomes better. 



IX. DCA IN FREQUENCY SPACE 

As illustrated for the momentum space, the DCA re- 
sults in a significant reduction of the problem complex- 
ity and it is complementary to the finite size lattice ap- 
proach. In analogy to the momentum space, one could 
consider dividing the one dimensional Matsubara fre- 
quency space into a number of coarse-graining subcells. 
For both fermions and bosons, each cell should include an 
odd number of frequencies in order for the frequencies in 
the centers of these cells to preserve Fermionic or Bosonic 
properties. Fig. 23 represents how the frequency space 
can be divided into coarse-graining subcells each com- 
prising a central fin frequency and a number of coarse- 
graining LVn lying around it. The central 0„ frequencies 
can be rewritten in the form of the original lattice with 



renormalized (3 shown as Pc iu the figure. Similar to 
the case of the momentum space, we make the following 
transformation for the Laue function 



o„i,o„2-i-r2„3 



(46) 



for the Matsubara frequencies of the vertex shown in 
Fig. 1 considering frequency dependent interactions in 
general (in condensed matter physics, most of the inter- 
actions are indeed simultaneous and thus frequency in- 
dependent). As a result, we may again coarse-grain the 
Green function over the subcell frequencies 

(5(K,f2„) = |^|]G(K,n„+w„). (47) 

According to Fig. 23, the full coarse- graining of the 
Green function amounts to /3c ^ which causes all the 
self-energy Feynman diagrams ordered higher than first 
(Hartree-Fock diagrams) to vanish and consequently we 
arrive at a fully static problem. 

Unfortunately, we can show that coarse-graining over 
Matsubara frequencies can lead to the violation of causal- 
ity and as a result, the DCA is not systematically imple- 
mentable for the Matsubara frequency quantities. The 
simplest example is the non-interacting Green function 
coarse-grained as follows [23] 

A 



(K, ^n) = jYl G'^°Hk, 0„ + Qn) = 



R 2-^ i 



1 



/3 4^ ifln + iWn - Cfe -h /X 



(48) 



The retarded Green function is derived by the substitu- 
tion fl-\-i'(] 



/? 4-^ O -h ir/ -I- ibJn - Cfe -I- /i ' 



(49) 



The theory of analytic functions of a complex variable 

tells us that in order for the Green function to remain re- 
tarded in the time space, G|^gj(K, ft) must not have any 
poles in the upper half plane of fl. In Eq. 49, one could 
readily create poles in the upper half plane for negative 
Hin which causes causality violation and consequently un- 
physical results. 
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FIG. 23. Dividing the fermion (a) and boson (b) Matsub- 
ara frequency spaces into subcolls. Tiio central fin froquoncios 
can be written in the form of the original lattice with renor- 
malized f3 denoted as /3c (/3c = 0/3 in this case). 



Solve Cluster 
Problem 



G(a) = 



1 



Q. + CO- ek+ H- T.'^iQ.) 



FIG. 24. The algorithm for the cluster solvers such as the 
NCA in which at each step before using the cluster solver 

the bare Green function is calculated by excluding the cluster 
self-energy. This is the only step through the algorithm where 
problems with causality might occur. 

The way around this difficulty is to use real frequencies 
even at finite temperatures. By invoking real frequencies 
we no longer coarse-grain over imaginary values. There- 
fore, the retarded Green function will never acquire poles 
in the upper half plane and remains causal as shown be- 
low 



Q + uj + ir] -\ — Efc + M — Sret(K, r2) ' 



(50) 



with N^^ the total number of frequencies and A^^ the num- 
ber of those we coarse-graining over in each cell. In an 
article by Hcttlcr et al. [4], a formal proof of causality 
is given (based on a geometrical argument) for coarse- 
graining in the momentum space. This proof can be 
straightforwardly applied to the real frequency space as 
well and it extends the application of real frequency DCA 
not only to the pcrturbative cluster solvers such as the 
FLEX but also techniques like the NCA. 

Lastly, similar to the momentum space, care must be 
taken when choosing the the size of the frequency coarse- 
graining cells. One must make sure that the cells are not 
larger than some characteristic energy scale (e.g. the 
Kondo temperature Tk) as the coarse-graining would 
then suppress the relevant physics. 



X. CONCLUSIONS AND OUTLOOK 

We introduce and examine the DCA in detail by em- 
ploying it with the FLEX to study the half filled two 
dimensional Hubbard model. The FLEX is not as pre- 
cise as nearly exact techniques such as quantum Monte 
Carlo in describing the Hubbard model at strong inter- 
action regime. However, it is capable of illustrating the 
utility of the DCA, including the complementarity and 
convergence of the DCA compared to finite size lattice 
approaches. The DCA and finite size calculations (with 
periodic boundary conditions) both converge with cor- 
rections 0{X/L^); however, in our example the coefficient 
^DCA was smaller and of opposite sign than \fs, indi- 
cating that the DCA converges more quickly and from 
a complementary direction. This complementarity was 
also seen in other quantities such as the pseudogap in 
the density of state and the non-Fermi liquid behaviour 
that the DCA (finite size) calculation systematically un- 
der (over) estimates. 

We also provide a detailed microscopic definition of the 
DCA by inspecting the the error generated by coarse- 
graining the Green functions in the compact and non- 
compact contributions to the thermodynamic potential. 
We conclude that due to the large magnitude of error 
that it generates, coarse-graining the Green function in 
non-compact part should be avoided and only the com- 
pact contribution should undergo coarse-graining. It also 
appears that coarse-graining the Green functions over the 
Matsubara frequencies can and will lead to the violation 
of causality and therefore is pathological. Nevertheless, 
one can coarse- grain the Green function over real frequen- 
cies and preserve the causality not only for the FLEX but 
also cluster solvers such as NCA in which the cluster con- 
tribution to the coarse-grained dressed Green function is 
excluded before being inserted in the cluster solver. 

The outlook for the FLEX-DCA approach is promising. 
Although the FLEX fails to accurately describe short- 
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ranged physics such as moment formation (and related 
phenomena hke the Mott gap) , it does a good job describ- 
ing long-ranged physics associated with spin and charge 
fluctuations. On the other hand, numerically exact cal- 
culations such as QMC, are too expensive to perform 
for large clusters, and are thus restricted to the study of 
short-length scales. However, since the DCA gives us a 
way of parsing the problem into different length scales, 
it may be used to combine the short-length scale infor- 
mation from the QMC with the long length scale infor- 
mation from the FLEX. This may be accomplished, by 
embedding a QMC cluster, of size L, into a much larger 
FLEX cluster of size L' ^ L, which is itself embedded 
in a mean-field. As we have shown here, this approach 
should be implemented by approximating the generating 
functional $ w $QMc(i) - '^flex{L) + <^flex{L'). 
Work along these lines is in progress. 
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